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Polariton condensates have proved to be model systems to investigate topological defects, as they 
allow for direct and non-destructive imaging of the condensate complex order parameter. The fun- 
damental topological excitations of such systems are quantized vortices. In specific configurations, 
further ordering can bring the formation of vortex lattices. In this work we demonstrate the sponta- 
neous formation of ordered vortical states, consisting in geometrically self-arranged vortex-antivortex 
pairs. A mean-field generalized Gross-Pitaevskii model reproduces and supports the physics of the 
observed phenomenology. 

PACS numbers: 03.65.Wj, 67.10. Ba, 71.36.+C, 42.50.Gy 



Quantized vortices are fundamental and ubiquitous en- 
tities across physics, playing a central role in mechanisms 
ranging from galaxy formation to phase conformation in 
microscopic quantum systems. They represent topolog- 
ical excitations of quantum degenerate Bose gases, such 
as Bose-Einstein condensates (BEC), superfiuids and su- 
perconductors pp. Throughout the last decades, these 
systems have offered the unprecedented opportunity of 
studying such topological defects and their phenomenol- 
ogy in a direct and controlled way [T]. Under peculiar 
conditions, quantized vortices have the unique property 
of arranging themselves in geometrically ordered struc- 
tures, such as the Abrikosov lattices [HE]. Vortex lat- 
tices were first observed in type-II superconductors under 
magnetic fields [4 a , then in both superfiuids [5] and atom 
BEC [MS], by setting the system into rotation, and also 
in optical nonlinear systems [101 E] ■ Ultimately, in the 
limit of high vortex density, these lattices are predicted 
to undergo a quantum phase transition to strongly cor- 
related states, similar to quantum Hall states, that still 
represent an open experimental challenge [lj. 

Recently, exciton-polaritons have established them- 
selves as a model two-dimensional Bose-gas [12]. Such 
quasi-particles arise as the eigenmodes of the strong- 
coupling regime between light and matter, which was 
demonstrated in planar semiconductor microcavities [T3] ■ 
The polariton system being dissipative in nature, due to 
the finite lifetime of the quasi-particles, its phenomenol- 
ogy is intrinsically and strongly out-of-equilibrium. 
Thus, continuous optical pumping is required to replen- 
ish the polariton population. Shaping of the excitation 
conditions allows manipulation of the condensate phe- 
nomenology in a simple way, unveiling striking physical 
effects. Moreover, thanks to the mixed light-matter com- 
ponents of polaritons, the emitted photons inherit all the 
properties of the quantum fluid that can, thus, be fully 



characterized by optical measurement of the extracavity 
field [H]. 

After the recent demonstration of polariton conden- 
sation [15 and superfluidity [THl [T7] in semiconductor 
microcavities, large efforts have been devoted to the 
study of quantum turbulence and vorticity, under dif- 
ferent excitation conditions, in such out-of-cquilibrium 
quantum fluids |18H2(Jj In particular, quantized vortices 
were demonstrated to spontaneously occur in the sys- 
tem as topological defects pinned by the disorder poten- 
tial |21j , which naturally results from the sample growth 
process. In analogy with the striking demonstrations 
of collective vortex phenomena in atom BEC, a num- 
ber of proposals and theoretical investigations considered 
means for the formation of vortex lattices in polariton 
condensates, both in the scalar [2"2~H2"5] and spinor con- 
densate case [231 121] • Methods for detecting the rotating 
lattice have also been proposed for the experimental ob- 
servation [2S]. Recent demonstration of vortex lattices 
were provided by exploiting periodic metallic structures 
that impose the lattice unity cell geometry [27] and by 
engineering the interference of multiple independent con- 
densates [2H] j bringing the topic to the forefront of po- 
lariton research. 

In this letter, we demonstrate the spontaneous occur- 
rence of polariton condensed states composed of geomet- 
rically self-arranged vortex-antivortex pairs. Each vortex 
having its counter-rotating partner, the overall state car- 
ries no orbital angular momentum. We perform interfer- 
ometric measurements allowing the determination of the 
amplitude and the phase of the condensate order param- 
eter, directly imaging the vortical entities. Theoretically, 
a mean-field approach reproduces the vortex-antivortex 
lattice formation. 

The sample is the same CdTe planar semiconductor 
microcavity, featuring 26 meV of Rabi splitting, used 
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FIG. 1: (a) Scheme of the setup used to engineer the ex- 
citation. A typical shaped intensity profile is shown in (b), 
averaged over many disorder realizations. Intensity profile for 
the pump beam used in the theoretical calculations (c). The 
red line depicts the cross-section intensity profile of the beam. 



in our previous experiments jT3J HH HH]- The sample 
is kept in a He-cooled cryostat at approximately 4 K. As 
mentioned, a key feature for the stabilization of vortex- 
antivortex lattices lies in the careful shaping of the in- 
tensity profile of the pump laser. The principle is to 
create a ring-shaped intensity profile, similar to a Bessel 
beam, with a large central spot of the same intensity 
as the outer ring. The reduced intensity valley between 
the central spot and the outer ring acts as the guide for 
the stabilization of vortex-antivortex pairs. Simply us- 
ing a Gaussian shaped pump would not be enough, since 
vortices are known to be dynamically unstable in such 
a configuration (3D]. The shaping of the beam is done 
by employing a single conical lens - axicon - in combi- 
nation with two lenses (LI, L2) and a 0.5 N.A. micro- 
scope objective (MO), according to the scheme depicted 
in Fig.[TJa). The initial top-hat laser beam is transformed 
into a Bessel beam by using the axicon. In the transfor- 
mation, intermediate intensity profiles are generated, as 
the one used in this work [Figjljb)], which is then im- 
aged by lenses LI and L2 in the real space of the sample, 
through the MO. Fine tuning of the excitation is allowed 
by a system of relative micromctric positioning between 
the axicon and the lenses. 

We excite the system above condensation threshold 
with the intensity-shaped pump beam at a power of 
250/iW/cm 2 The photoluminescence (PL) coming from 
the sample is collected by the same MO and is then sent 
to a modified Mach-Zehnder interferometer, which takes 
the condensate image and interferes it with a 4x magni- 
fied replica the central spot (marked by the dashed rect- 
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FIG. 2: (a) Polariton condensate density at a position where 
a four-vortex lattice is observed, (b) Magnified central part 
of the condensate [marked by the dashed rectangle in panel 
(a)] used to overlap with the full condensate to obtain the 
interference pattern, shown in (c). The corresponding phase 
structure is plotted in (d). Red (green) circles mark the po- 
sition of vortex (antivortex) . 



angle in Fig. [2ja)). This allows to take a small part of 
the condensate, acting as phase reference, overlap it with 
the whole condensate image and extract its phase struc- 
ture, as performed in Ref. [29 . From the interference 
pattern, the phase can then be extracted by digital off- 
axis holographic techniques [31] • Despite the presence of 
a disorder potential [32] , some regions of shallow disorder 
are present in the CdTe sample that allow the formation 
of vortex lattices. 

To explain the whole process, we first consider a set 
of data showing a four-vortex lattice, composed of 2 
vortex-antivortex pairs. The condensate density image 
(Fig. |2ja)) is interfered with the enlarged version of the 
central spot of the condensate (Fig. [2jb)). The resulting 
pattern is shown in Fig. [2jc) together with the corre- 
sponding phase structure reported in Fig. [2jd). In the 
figure, two vortex-antivortex pairs can easily be identi- 
fied, marked by the red (vortex, clockwise phase winding) 
and green (antivortex, counterclockwise phase winding) 
circles at the position of the phase dislocation. Observ- 
ing the real space phase map, when circumventing the 
center of the condensate cloud, one finds alternately vor- 
tex/antivortex entities. Let us point out that some spu- 
rious phase singularities (with no circular marker on top) 
appear in a region of negligible density, outside the con- 
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FIG. 3: (a) Polariton condensate density at a position where 
a six-vortex lattice is observed, (b) Profile taken around the 
central spot of the condensate that highlights the relative min- 
ima, identifying the core of the six vortices (each vertical line 
marks the minimum position), (c) Interference pattern from 
which the phase structure is extracted, as shown in (d). Red 
(green) circles mark the position of vortex (antivortex) . 



dcnscd polariton gas. They are most probably resulting 
from overlap of the decaying outer tails of the condensate 
with the edges of the flat phase reference. 

It is possible to create vortex lattices of higher order, 
where more vortex-antivortex pairs are able to accom- 
modate in the same alternate ordered way. By slightly 
increasing the size of the pump spot (in the one micron 
range) while keeping its intensity profile as in Fig. [TJb), 
we could find other positions on the sample where the 
disorder allowed the formation of the six-vortex lattice. 
The results are summarized in Fig. [3j The condensate 
density, featuring six clear density minima is shown in 
Fig. [3ja).The vortex cores are evidenced also by taking a 
circular profile of the density [33], shown in Fig. [3jb). 
The vertical dashed lines identify the position of six 
clear relative density minima, vortex cores, along the az- 
imuthal profile. By analysis of the orientation of the fork- 
like dislocations in the interference pattern [Fig.[3jc)] and 
corresponding phase structure [Fig. |3jd)], one can iden- 
tify even better the position of the phase dislocations 
and the winding sign of each vortex. The red and green 
color circles, as in the case of Fig. [2] mark the alternating 
vortex-antivortex lattice. Despite the presence of the un- 
derlying disorder potential felt by the polaritons, the six- 
lattice nicely matches an hexagonal structure. Indeed, 
the disorder accounts for the distortions with respect to 





FIG. 4: Theoretically calculated polariton condensate den- 
sity (a) and phase structure (b) that nicely matches the ex- 
perimental results of Fig. [3] 



an exact symmetric, theoretically expected, Abrikosov 
vortex-lattice observed in other systems [SHIT]. 

Polariton condensation under non-resonant pumping 
is a highly non-equilibrium process. Being coherent, the 
condensate itself can be treated with a Gross-Pitaevskii 
type equation coupled to a classical rate equation for the 
dynamics of an exciton reservoir |34j . which appropri- 
ately models the incoherent excitations generated by non- 
resonant pumping that feed and interact with the conden- 
sate. Both spectral and spatial features coming from ex- 
periments have been successfully reproduced within this 
theoretical framework [35rf5?] and also adapted to take 
into account the dynamic properties of polariton con- 
densates [55rf4"0] and of spinor condensates [H][52]. The 
equations read: 



ih 



dt 

dn(r, t) 
dt 



ih 



E LP + V{v) + -{R R n{v,t)- lc ) 



(1) 

= -(lB + Rr\1>(t, t)\ 2 ) n(r, t) + P(r) (2) 



where ^(r, i) is the polariton condensate in the mean- 
field representation and n(r,t) is the intensity distribu- 
tion of the incoherent excitonic reservoir, directly in- 
jected by the non-resonant pump. For the purpose of 
this work the polarization degree of freedom of polari- 
tons is neglected. Elp is the polariton kinetic energy 
operator, accounting for the non-parabolicity of the dis- 
persion of lower branch polaritons. V(r) represents an 
effective potential, given by the mean-field shift induced 
by polariton-polariton interactions (g), the interaction 
of polaritons with the reservoir (go) and an additional 
pump induced shift (G) [3J]; V(r) = hg\tp(r,t)\ 2 + 
hg R n(r,t) + hGP(r), where g, g R and G are constants. 
P(r) represents the spatial pump distribution, as in 
Fig.[ljc). The intra-condensate interactions between po- 
laritons being present but small, the dominant effects 
determining the induced effective potential come from 
the effect of the exciton reservoir on condensed polari- 
tons [35] • 7c and jr represent the decay rates of con- 
densed polaritons and reservoir excitons, respectively. 
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Rr is the stimulated scattering rate of excitons from the 
reservoir to the condensate and thus represents the con- 
densation rate in the system. The system of two coupled 
equations [T] and [2] is numerically solved in time. The 
initial conditions for seeding the condensate consist of 
an initial random low-intensity white noise. Theoreti- 
cally, the system reaches the same steady vortex lattice 
pattern independent of the initial noise. The initial noise 
may however affect the time required for formation of the 
lattice. One can note that the vortex lattice pattern is 
not cylindrically symmetric (although it has a rotational 
symmetry of third order). Theoretically, the breaking 
of cylindrical symmetry could be caused by the initial 
noise, in which case the pattern would form with a dif- 
ferent overall orientation (but otherwise identical shape) 
from shot to shot. Our system is however not perfectly 
symmetrical to begin with; imperfections in the pump 
profile and disorder already break the system symmetry 
allowing the pattern to form with the same orientation 
from shot to shot and be observed in multi-shot averaged 
experiments [25] . 

Figure [4] shows the intensity and phase of the po- 
lariton condensate steady state, calculated with this 
method [43] under the custom intensity-shaped pump 
beam of Fig. [TJc) . A good agreement is found between 
theory and experiment (Fig. [4} both in the condensate 
density and the phase structure, which feature the three 
vortex- anti vortex pairs. The excitation laser, with its 
reduced intensity central ring, acts as a guide and as a 
trap for the vortices and favors the formation of spatially 
coherent states with regular lattice of vortex-antivortex 
pairs within the trapping potential due to its shape. 

The number of vortices composing the lattice is partly 
sensitive to the profile of the pump, in terms of its size 
and shape. The use of larger sized pumps allows for larger 
vortex lattices with more vortex-antivortex pairs. Fur- 
thermore, it is worth mentioning that the disorder plays 
an important role in the determination of the energy 
of the condensate, substantially influencing the poten- 
tial landscape felt by the particles, so that only at a few 
specific sample positions could we obtain the engineered 
pump-induced trapping of the vortices and observe the 
lattice formation. At such positions, slight adjustment 
of the pump beam alignment allows fine tuning of the 
symmetry and stability of the lattice itself. 

The theoretical results presented in Fig. [4] show that a 
steady-state has formed in the system for the chosen form 
of excitation. However, it should be appreciated that 
steady-state solutions may not exist for all pump profiles 
and intensities. For a smaller sized-pump we have also 
obtained periodic solutions of a rotating vortex lattice, 
where the number of vortices and antivortices need not 
be an even number. The experimental observation of 
such states is a challenging task, either requiring single- 
shot time-resolved interferometry or by means of more 
advanced interferometric techniques [26 . 



In this work we have reported on the formation and 
stabilization of vortex lattices, composed of vortex- 
antivortex pairs, in a polariton condensate. Ordered pat- 
terns are observed under proper intensity-shaped non- 
resonant optical pumping. The shaping is such to cre- 
ate a trap for vortices that are seen to arrange in a ge- 
ometrical vortex lattice. A delicate interplay between 
the excitation shape and the underlying disorder poten- 
tial pins the orientation of the spatial patterns, allow- 
ing their detection in time-integrated experiments. We 
make use of a mean-field approach in the form of a gen- 
eralized Gross-Pitaevskii model that is able to reproduce 
the lattice formation under the same kind of excitation 
conditions used in the experiments. The control over the 
lattice formation could prove to be useful in Bose gases 
when considering the few particle per vortex limit, which 
is expected to give rise to quantum phase transitions to 
highly correlated states [T] - a world yet to be explored. 
Recent long-lifetime polaritons achieved in GaAs cavi- 
ties [44] represent to this extent really promising systems. 
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